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RECONSTRUCTION OF A POLYNOMIAL FROM ITS RADON 

PROJECTIONS 

BORISLAV BOJANOV AND YUAN XU 


Abstract. A polynomial of degree n in two variables is shown to be uniquely 
determined by its Radon projections taken over [n/2] + 1 parallel lines in each 
of the (2[(n + l)/2] + 1) equidistant directions along the unit circle. 


1. Introduction 

Let / be a function defined on the unit disk B 2 on the plane. A Radon projection 
of / is the integral of / over a line segment inside B 2 . More precisely, for any 
given pair (9,t) of a real number t £ [—1,1] and any angle 9, let I(9,t) denote 
the line segment inside the unit disk B 2 , where the line passes through the point 
(t cos 0, i sin (7) and is perpendicular to the vector (cos 9, sin 9). Then 

(1.1) K e (f;t) := [ f(x,y)dxdy 

defines the Radon projection of / on the line segment 1(9, t). 

The Radon transform / i—> {lZg(f;t)} associates to / a family of univariate 
functions of t, parameterized by 9. The problem of reconstructing / from full 
or partial knowledge of {7 Zg(f;t)} has been studied by many authors. It plays 
an essential role in the computer tomography. It is well-known that the set of 
Radon projections {7 Ze(f;t) : 0 < t < 1,0 < 9 < 27t} determines / completely. 
Furthermore, it is known that if / has compact support in B 2 , then / is uniquely 
determined by any infinite set of Radon projections m, In practice, however, the 
data set is usually finite. Thus, the main problem is to obtain a good approximation 
to the function from a large collection of its Radon projections. See m for the 
background and detail discussions. Using a finite data set, one may determine a 
polynomial whose Radon projections agree with the given data. Such an approach 
has been considered in mm\m in connection with computer tomography. 

The present paper concerns with the problem of whether a polynomial of degree 
n in two variables can be uniquely determined from a set of (n+l)(n + 2)/2 distinct 
Radon projections. The solution depends on the arrangement of the lines on which 
the projections are taken. One solution was given early in ED, in which the Radon 
projections are taken over all possible line segments [x^, x^] joining x* and x^, where 
xq, xi,..., x n+ i are n+2 equally spaced points on the boundary of the unit disk B 2 . 
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In ^ (see also 0), a more general result was proved, in which xo, xi,...,x n +i are 
n + 2 distinct points on the boundary of a convex set and the result was extended 
in |7| to higher dimensions. Recently this question was considered in |T], in which 
the lines are taken in (n + 1) directions, with one line in the first direction, two 
lines in the second direction, etc. This set has the drawback of lacking symmetry. 

The main result of the present paper shows that a set of Radon projections taken 
over 2[n/2j + 1 parallel lines on each of the 2|_(n + 1)/2J + 1 equidistant directions 
can determine a polynomial of degree n uniquely. The set of the corresponding line 
segments possesses a rotation symmetry. If n = 2m is a fixed positive integer, then 
our method requires the Radon data in 2m + 1 directions, which are taken to be 
in equally spaced angles along the unit circle, and in each direction we take m + 1 
parallel lines, associated to t 0: ti,..., t m . We prove that for almost all choices of 
{tfc} a polynomial P n of degree n = 2m in two variables is uniquely determined by 
the set of its Radon projections over these line segments; more precisely, for any 
given /, P n is uniquely determined such that 

(1.2) TZ^(P n ;t k ) = TZ (/>j {f;t k ), 0 < j < 2m, 0 < k < m. 

The similar construction works for n = 2m—1. Several examples of the sets of points 
{t k } that define regular interpolation are given. Furthermore, the polynomial P n 
can be easily computed. Thus, the result appears to offer a simple way to recover 
a polynomial from its Radon data. 

The number of conditions o is equal to (n+l)(rc+2)/2, which is the dimension 
of the space of bivariate polynomials of total degree at most n. Hence, C3 can be 
interpreted as an interpolation problem by polynomials based on line integrals. As 
in the case of interpolation based on points, this bivariate problem is not always 
regular. The only general configurations in the literature that hold for every integer 
n are that of Hakopian 0 mentioned above and the non-symmetric configuration 
given recently in [Tj. The result of the present paper provides another family 
of examples. The proof is based on an observation that recovering polynomials 
from their Radon projections can be reduced, using a formula in in, to a family 
of univariate interpolation problems that uses certain special classes of algebraic 
polynomials. In the case of pQ, the strategy of solving the interpolation problem 
resembles the approach that uses the Bezout theorem for pointwise interpolation. 
The approach in the present paper follows the strategy of (2], which uses equally 
spaced points on circles for pointwise interpolation, that allows one to step beyond 
the limitation of the Bezout Theorem. 

The paper is organized as follows. In the following section we prove the existence 
and the uniqueness of the reconstructing problem. In Section 3 we show how to 
construct the polynomial and give an outline of the algorithm. 

2. The existence and the uniqueness of the solution 

2.1. Preliminary. Let B 2 = {( x,y ) : x 2 + y 2 < 1} denote the unit disk on the 
plane. Let 9 be the angle in the polar coordinates 

x = rcos9 1 y = rsm9, r > 0, 0 < 9 < 2n. 

A line l whose slope is — cot 9 is defined by the equation 
i(x, y) := x cos 9 + y sin 9 — t = 0, 
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where t is a real number. Clearly the line i passes through the point (t cos 9,1 sin 9) 
and is perpendicular to the vector £ = (cos 9, sin 9). Since, for a fixed t, the line 
corresponding to ( 9 , t) coincides with the line corresponding to (tt + 9, — t), we could 
assume 9 £ [0, tt). Alternatively, we could assume that 9 £ [0, 2tt) and t > 0. 

We will also use l to denote the set of points on the line £ and introduce the 
notation 

1(9, t) = £UB 2 , 0 < 9 < tt, -l<t<l, 

to denote the line segment of £ inside B 2 . The points on I(9,t ) can be represented 
as follows: 

x = t cos 9 — s sin 9, y = t sin 9 + s cos 9, 
for s £ [— Vl — t 2 , \Jl — t 2 ]. 

The Radon projection of a function / in the direction 9 with a parameter t £ 
[—1,1] is denoted by lZg(f;t), 


Kg(f;t) := [ f(x,y)dxdy 
Jue.t) 


'1(6,t) 

/Vi-t 2 

= / f(t cos9 — ssm9,tsm9 + s cos9)ds. 

J-VUU 2 

In the literature it is also called an X-ray. It is clear that 1Z$(f;t) = lZ n+ g(f; —t) 
so that we can assume, for the sake of definiteness, either O<0<7rorO<t<l. 


2.2. Polynomial bases. Let 11^ denote the space of polynomials of total degree n 
in two variables, which has dimension dim 11^ = (n + 1 )(n + 2)/2. If P £ 11^ then 

k 

x = (x,y). 

k— 0 j —0 

Let Uk denote the Chebyshev polynomial of the second kind, 

sin(fc + 1)9 


p ( x ) = 


Uk(x) = 


sin 9 


x = cos 9. 


For £ = (cos 9, sin 9) and x = (x,y), the ridge polynomial Uk(9; ■) is defined by 

U k (9;x) := C4((x,£)) = U k (xcos9 + ysinfl). 

Clearly Uk is an element of 11^ and it is constant on every line that is perpendicular 
to £. The Radon projection of this function in any direction can be easily computed. 
This is a result due to Marr cn and it plays a central role in our discussion below. 

Lemma 2.1. For each t £ (—1,1), 0 < 9, </> < 27r, 

7l4U k (9; •); t) = -Vl-t 2 U k (t)U k (cos(<f> - 9)). 

The following useful relation follows easily from Marr’s formula m), 


( 2 . 1 ) 


- / U k (9;x)U k ((t>;x)dx = —— U k (cos(cj) - 9)). 
ttJ B 2 k +1 


Let V n denote the space of orthogonal polynomials of degree n on B 2 with respect 
to the unit weight function; that is, P £ V„ if P is of degree n and 


[ P(x)Q(x)dx = 0, 
J B 2 


for all Q £ II^_ 1 . 
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For special choices of 9 and <f >, the equation EH becomes an orthogonal relation. 
Indeed, for k £ N, let 

j 7T 

fj,k — (cos sin ^ ? 0 + 3 — 

For fixed k, the points cos 6j tk , 1 < j < fc, are zeros of U k ■ The ridge polynomials 
£4(0; x) have the remarkable orthogonal property that Uk{6 ; •) £ 14. Since the 
dimension of 14 is k + 1, it follows from ED that the set 

Pfc := {U k {9 jtk ;x) : 0 < j < k} 

is a basis for V k . In particular, this shows that the set {Pfc : 0 < k < n} is a basis 
for Il£. Together with Lemma. 12. II this proves the following result: 

Lemma 2.2. Every polynomial P n £ 11^ can be written uniquely as 

n k 

(2-2) W=££ 

k =0 j =0 

Furthermore, for each </> and t, 

n o k 

(2.3) TZ^(P n ;t) = \/l - t 2 £ —-14(1) £cj, fc I/ fc (cos(0 - gj, fc )). 

/c—0 j=0 

There are several other explicit orthogonal bases for V n ; see, for example, 21. If 
{Q k : 0 < k < n} is an orthogonal basis of V n , then it was shown in m that 

k 

(2.4) U k (9j }k ;x) = —^—'^2Qf(x,y)Qf(cos9 jtk ,sm9 jik ). 

+ 1 ;=o 

One explicit basis of V n , denoted by <3j,i( x ), is given in terms of polar coordinates 
as follows (cf. 0), 

Qj,i(x,y) = h n - 2 j,iPj°' n ~ 2:,) (2r 2 - 1 )r n ~ 2j cos (n - 2 j)8, 0 < 2 j < n, 

Qj^(x,y) = h n - 2 j, 2 Pj°’ n 2j \2,r 2 - 1 )r n ^ 2j sin(n - 2 j)0, 0 < 2j < n - 1, 

where P denotes the usual Jacobi polynomial, ho,i = l/(n+ 1)) and hj t i = 
hj 2 = 1/(2 n + 2) for 1 < 2j < n — 1. Using this basis and restricting x to the 
boundary of B 2 , the equation 12.411 becomes 

k 

U 2 k (cos(</> — 9j 2 k)) = 1 + 2 (cos 219j 2 k cos 2l(f> + sin2l9j t 2 k sin2Z</>), 

1=1 

(2-5) * 

U 2 fc-i(cos(0 - 9 jt 2k- 1)) = 2 £ (cos(2 1 - l)0j, 2 fc-i cos(2Z - 1)</ 

i=i 

+ sin(2Z — l) 0 j, 2 fc_i sin(2Z — 1 )</). 


The above equations also follow from the elementary trigonometric identities 

sin(2fc + 1)0 £ sin2fc0 

1 + 2£,cos2j0 = -t— - and 2 £ cos(2j — \)9 = 


J=i 


l=i 


sin 0 
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with 6 replaced by <fi — d^ 2 k and </> — Oj^k- 1 , respectively. These relations allow us 
to rewrite m in a form that is easier to work with. In the following we use L^J 
to denote the integer part of x. 


Lemma 2.3. Let P n be given as in ( 12.211 . Then 

Ln/2J 

( 2 . 6 ) 


~7T = W 


= E ^j(t) 


1=0 


ao. 2 i + 2 E ( a i,2i cos 2 /</> + bj, 2 i sin 2 j<jj) 


i=l 


L(n+i)/aj 
z=1 


2 E( Q i, 2 Z-i cos(2j - l)</> + 6 j, 2 Z-i sin(2j - 1 ) 0 ) 
i=i 


m which the coefficients Cj t k in (EH and aj t k, bj'k are related by 
(2.7) Cj, 2 i = —« 0 , 2 Z + ^2 (afc ,2 1 cos 2 kOj^i + h, 2 i sin 2k0j^i) 


k =1 


for 0<j<2l,0<l< [n/2\, and 
i 

(2.8) Cj,2i-1 = E («fe,2Z—l cos(2fc - l)6»j,2Z-i + 6fc,2i-i sin(2fc - 1)6 » Jj2 z—1) 


fc=i 


/or 0 < j < 21 - 1, 1 < l < [(n + 1)/2J. 


Proof. We only prove the case of n = 2 to. Comparing (T2~3l) and (ESj) shows that 


1 l ; 

— E c i,2Z^2z(cos((/ - 0j t 2i)) = ~ a o,21 + E( a k 2i cos2 // > + & J,2Z sin 2 j 0 ), 


2Z + 1 


i=o 


i=i 


2Z-1 


yj E C J,21-1^21-1 (cos(0 — 


21 


j =o 


= E^J- 2 '- 1 C0S ( 2 / - 1 ) < ^ > + fo J,2Z-l sin ( 2 J - 1 )</>)■ 
j=l 

Using the equations in 12.511 we can rewrite the left hand side in the form of the 
right hand side, which leads to the equations: 


21 


21 


2 x 2 x 

aj,2i = 21 i i / . c fc,2Z cos 2jdk,2i, bj t2 i = . } y Cfc^z sin 2 j0k,2i 

' 1 fc =0 + 1 fc =0 

for 1 < j < l and 2 ao ,22 satisfies the above formula with j = 0, and 

, 21 -1 


j,2Z-l = y E C MZ-1 cos (2 j — l)9k,2l-U 


fc=0 
2Z-1 


bj, 2 i-i = y E sin(2/ — l) 0 fe, 2 Z— l 


fc =0 


for 1 < j < l. For fixed Z, the above linear relations can be written in the matrix 
form. After a proper normalization, the coefficient matrix turns out to be an 
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orthogonal matrix. For example, for the even indices, we have 



r 1 

1 

1 



0-0,21 


72 

72 COS 

72 cos 


co ,21 


a 1,2 1 

y* 

1 

0 


72 sin 


Cl ,21 

1 

01,21 

V2 

V2Z + 1 

72 

72 cos 2 ••• 

72 cos 2 


C2l,2l_ 

y/2l + l 

0-1,21 

V2 


0 

72 sin 2Z 

72 sin 21 ffjzpi 



°l,2l 

. V2 . 


The (l + 1) x (l + 1) matrix in the left hand side, including the factor l/y/2l + 1, 
is orthogonal. Hence, the linear system of equations can be easily reversed, which 
leads to the stated relations E3 and (El). The formulas can also be verified 
directly by inserting E 3 and E 3 into the equations of and b 3) k and using 
the well-known trigonometric identities such as 


1 

m + 1 


£cOS 2j0k,m 

k =0 


1 


m + 1 


^ ^ sin 2 j6k,m 
k—0 


m + 1 

o 


e 2ij6k - m = 6 , 


'o ,j 


for 0 < j < to, where Sqj = 1 if j 7 ^ 0 and Sqj = 0 otherwise, and the elementary 
trigonometric identities such as 2 cos 9 cos (j> = cos(f? + (/>)+ cos (9 — (f). □ 


To reconstruct the polynomial, we will determine the coefficients and byfc and 
use them in E 3 > and E 3 ) to determine Cj,k in ESI. The following representation 
is useful. 


Lemma 2.4. Let n = 2m or n = 2m — 1. Let P n be given as in 12.21 . Then 


(2.9) 

where 




and for 1 < j < m, 


n 

A 0 (t) + Y [Ajit) cos j</> + Bj(t) sin j<j >\, 
3 =1 

L n /2J 

Ao(t) = Y ao,2iU2i(t), 

1=0 


( 2 . 10 ) 

( 2 . 11 ) 


Proof. 


[n/2J ln/2J 

A2j(t) = 2 Yj a j,2lU2l{t), B2j(t) = 2 Y^ bj t 2lU2l(t), 

1=3 1=3 

L n /2J h/2J 

A2j-i{t) = 2 aj i 2i-iU2i-x{t), B2j-\{t) = 2 "Y^ ^3,21- 1U21-1 ( i ). 
1=3 1=3 

This comes from changing the order of summations in itroli . □ 


2.3. Existence and uniqueness of the solution. We now state the Radon pro¬ 
jections from which the polynomial P n can be uniquely recovered. These are given 

by 

( 2 . 12 ) ^, m (P„;4), 0 < fc < |_n/ 2 j, 0 < j < 2m, m=[(n+ 1 )/ 2 J, 

which consists of total (|_n/ 2 j + l)(( 2 |_(n + 1 )/ 2 J + 1 ) = (n + l)(n + 2)/2 many pro¬ 
jections, the same as the dimension of n^. The angles are chosen to be equidistant, 
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(2.13) 6 m := | ^ : 0 < j < 2m|. 

The reason that we choose the equidistant angles lies in the lemma below, which 
plays a key role in our study. Such a lemma appears first in |2j and has been used 
in PJ and ITTj . 


Lemma 2.5. For <fi £ 0 m and P n £ 11^ with n = 2m or n = 2to — 1, 
T^<j>{P n \ t) 




= A 0 (t) 


+ ^ [(Aj(*) + cosjcj) + (Bj(t ) - B 2m _j + i(t)) sin j(j>] , 

j=i 


where we assume that A 2m = I? 2m = 0 if n = 2to — 1. 


Proof. Using the expression in the previous lemma, the proof follows from the fact 
that 

cos(2 m — j + 1 )</> = cos j(f> and sin(2m — j + 1)4> = — sin j<f 
for cj) £ 0 m . □ 


For the expression in this lemma, the variable t is in (—1,1). We can also state 
the result for </> £ 0 m = {(2 j + 1)7t/(2to+1) : 0 < j < 2to}, for which the sign before 
A‘ 2 m -j+i and will be reversed. It is easy to see that 0 m = 0 m +7r modulus 

27 t . The two expressions are consistent with the fact that TZg(f\t) = IZ^+t »(/; —t). 

We are now ready to prove our main result, which shows that the polynomial 
P n can be uniquely determined by the data (E~m . We need the following function 
classes: For n = 2m define 

Xj(t) = {U2m(t), U2m—2{f ); • • ■ , U2j(t), U2m-l(t), U2 m s{t), . . . , U2m-2j+l (t)} 

for 1 < j < m — 1. We also consider Xj as a column vector in K. m+1 and regard 

= [Xj(t o)i ... 5 Xj(tm)\ 

as an (to + 1) x (to + 1) matrix. 

A set of points t = {to,ti, ■ ■ ■ ,t m } in (—1,1) is said to be asymmetric if t and 
—t do not both belong to t. 

Theorem 2.6. Let t := {to,Uj • • - Am} be a set of asymmetric distinct points in 
(— 1 , 1 ) such that the matrices 5i(t), S 2 (t),..., S m _i(t) are all nonsingular, then 
the polynomial P £ 11^ is uniquely determined by the set 12.1211 of its Radon 
projections. 

Proof. In order to prove that there is a unique solution, it is sufficient to show that 
if 72.0 i m (P; tk ) = 0 for 0 < j < 2m and 0 < k < m, then P(x) = 0. 

Setting its left hand side zero, the expression in Lemma 1‘2 , 51 shows that 

m 

A 0 (t k ) + "^2 [(Aj(tk) + A 2m -j+i(tk)) cos j(j) + ( Bj(tk ) - B 2m -j+i(tk)) sin j<f\ = 0 
3 =1 
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for 0 < k < m and </> £ @ m . The left hand side is a trigonometric polynomial of 
degree m and it vanishes on 2m + 1 points, hence, it follows from the uniqueness 
of the trigonometric interpolation that the coefficients have to be zero; that is, 

Ao(tk) = 0, Aj(tk ) + A-2m—j+l(tk) = 0, 

Bj(tk) - B 2m -j+i(t k ) = 0, 1 <j<m 

for 0 < k < m. Since Aj and A 2m -j +1 have different parity, the coefficients in 
Aj(t) + A 2 m -j+i{t ) will not combine and there are exactly m + 1 terms in this 
polynomial, written as a linear combination of the Chebyshev polynomials {Uk}- 
Let Y 0 (t) = {U 2 m{t), U 2 m- 2 (t), ■■■, U 2 (t ), U 0 (t)} and let 

Y 2 j(t) = {U 2m {t ), U 2rn - 2 (t), . . . , U 2 j{t ), U2m-3(t), ■ ■ ■ , U2m-2j+l{t)}, 

Y 2 j-l(t) = { U2m—1 (t), U 2 m-3(t), ■ • • , U 2 j+l(t), U 2 m{t), U 2m - 2 (t), ■ ■ ■ , U 2m - 2 j(t)}. 

Furthermore, let Sy^ denote the matrix 2y^. = [Y)(to), Yj (ti),..., Yj(t m )}. The 
coefficients of the linear systems of equations in 12.141 are 2y 0 , 2y 2j . for 1 < j < 
m/2, and Sy 2j _j for 1 < j < (m + l)/2. Hence, in order to prove that 12.141) 
implies H 0 (t) = 0, Aj(t) = Bj(t) = 0 for 1 < j < m, we need to show that these 
matrices are all invertible. However, it is easy to see that we have the relation 
Sy 2(m _..)_i = Sy 2j , which shows, in particular, that Sy m = Hy m _ x . Thus, we can 
deal with Sy 2i for 0 < j < m — 1. Furthermore, Yq contains only even polynomials 
U 2 i, 0 < l < nn. Using the notation s = t 2 , it follows that U 2 i{\fs) is a polynomial 
of degree l, so that the matrix Sy 0 is always invertible as t is a set of asymmetric 
distinct points. Thus, since 2 y 2j = 2j(t). our assumption on t implies that all 
matrices are invertible. □ 

A similar theorem holds in the case of n = 2m — 1. First we need to define for 
n = 2 in — 1, 

Xj (t) = {U2m- 2 (t), U2m-4:(t), ■ ■ ■ , U2j(t), U2m-l(t), U 2rn s(t), . . . , U 2m -2j+l(t)} 
for 1 < j < m — 1 and 

X m (t) = {U2m-l(t), U 2 m-3(t), ■ ■ ■ > U 2 (t), Ui(t)}. 

We will keep the notion 2j (t) for the matrix built upon from the columns of Xj, 
which is an m x m matrix. 

Theorem 2.7. Let t = {to, ti,..., t m -i} be a set of asymmetric distinct points in 
(—1,0) U (0,1) such that the matrices Si(t),S 2 (t),... ,2 m (t) are all nonsingular, 
then the polynomial P £ n 2m -l * s uniquely determined by the set (I2T21) of its 
Radon projections. 

That the points to, ti,..., t m are all in (—1, 0) U (0,1) means that none of the 
points can be zero. In fact, since U 2 k-i is an odd polynomial, 2 m (t) contains 
only odd polynomials, which vanish at the origin. Otherwise, the proof of this 
theorem is similar to that of the previous theorem. In the case of n = 2m, the 
set Xq = {U 2 m, U 2m - 2 , ■ ■ ■, Uq} does not appear in the conditions of the theorem 
since the interpolation at distinct points in [0,1) by Xq is always regular. Such a 
reduction of conditions does not appears to happen for n = 2m — 1. 

For each fixed j, the determinant of 2 ? (t) can be consider as a polynomial 
function in t = (to, t ±,..., t m ), so that det2j(t) = 0 defines a hypersurface of 
dimension m. Hence, the determinant is not zero for almost all choices of t £ R m+1 . 
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Evidently, the same holds true for m + 1 matrices. Hence, we have the following 
corollary: 

Corollary 2.8. For n = 2m or 2m — 1 and for almost all choices of distinct 
points to,ti ,..., t m in (— 1 , 1 ) for n = 2m or in (— 1 , 0 ) U ( 0 , 1 ) for n = 2m — 1 , 
the polynomial P n £ n 2 is uniquely determined by the set lt 2 J 2 l) of its Radon 
projections. 

2.4. Choices of tk . We have shown that almost all choices of t = {to, ti,..., t m } 
will lead to the unique solution of recovering the polynomial. One may ask the 
question if any choice of t C (—1, l ) m+1 will work. The answer, however, is negative 
as the following example shows. 

Example: m = 2. By Theorem 12.61 we only have to choose to,ti,t- 2 in (0,1) such 
that the set X 2 = {U 2 , U 3 , t/ 4 } has nonsingular determinant 3 2 (t). Recall that 

U 2 (t) = 4t 2 - 1, U 3 (t) = 8 f 3 - 4i, U 4 (t) = 16f 4 - 12t 2 + 1. 

We can compute the determinant 32 (t) explicitly and the result is 

det S 2 (t) =32 (ti - tj) 

l<i<j<3 

x [ 8 + 4tit 2 t 3 (ti + t 2 + t 3 ) + ( 2 1\ - l)(2tl - l)(2tg - 1)] . 

Since ti £ (0,1), the first two terms in the square bracket are positive, while the 
third term could be negative. This shows, in particular, that det S 2 (t) 0 if one of 

the ti is v / 2/2 or if two of ti,t 2 ,t 3 are less than V2/2 and the other one is greater 
than v / 2/2. However, detS 2 (t) can be zero for some choices of to. t\. f 2 . 

On the positive side, we give two results that provides sets of points that will 
ensure the uniqueness of the reconstruction. The first result uses the following 
theorem due to Obrechkoff Ill- 

Lemma 2.9. Let dp be a nonnegative weight function on an interval and let Pq, Pi, 
P 2 , ... , be the orthogonal polynomials with respect to dp. Denote by a n the largest 
zero of P n . Then the number of zeros of the polynomial oqPq + diP 2 + ■ ■ ■ + a n P n , 
where do, a\,..., a n are real numbers, in the interval (a n , +00) is at most the num¬ 
ber of sign changes in the sequence of the coefficients ao,a±,... ,a n . 

The Chebyshev polynomials Uo,Ui,U 2 , ■ ■ ■ are orthogonal with respect to the 
weight function -\/l — x 2 on [— 1 , 1 ], so that the above lemma implies the following: 

Proposition 2 . 10 . Let t = {to, ti,..., t m } be a set of numbers that satisfies 

(2.15) cos n <t 0 <ti < ... <t m <1. 

2m + 1 

Then the matrices 5i(t),..., S m _i(t) are all nonsingular. 

Proof. If one of the matrices, say 3 ? - (t), were singular, there would be a nonzero 
polynomial P £ span Xj that vanishes on t. This means that the number of zeros 
of P in (a n , 00 ) would be m + 1. However, each set 3 j has cardinality m + 1 so 
that the number of sign changes in the sequence of the coefficients of P is at most 
m. This is a contradiction to the conclusion of the lemma. □ 
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The set of points given in this proposition ensures the uniqueness of determining 
a polynomial by the set of its Radon projections. The condition 112.1511 implies that 
all points are clustered toward the one end of the interval [—1,1]. This appears to 
be neither practical nor a good choice for computation. In fact, the numerical test 
indicates that some of the matrices Hj- tend to have very large condition numbers. 

Our second positive result is more interesting. Here the points {tfc} are based 
on the zeros of the Chebyshev polynomials [/ 2m . It is well-known that these zeros 
are given by 

J7T 

%-,2m := C0S 2^r[i j = 1,2,...,2m. 

We state the result for n = 2m first. 

Theorem 2.11. Let to be any point in (—1,1) such that f7 2m (to) 0- Let 

tj = V 2 j, 2 m for j = 1,2 ,...,m. Then the matrices S 1 (t),S 2 (t),... ,S m _i(t) in 
Theorem 2.6 are all nonsingular. Consequently, the polynomial P £ is uniquely 
determined by the set itm of its Radon projection. 

Proof. It is easy to see that the set t is asymmetric. Hence, according to Theorem 
2.6, we only need to show that Ej(t) is invertible for each j = 1, 2,..., m— 1. Using 
the explicit expression of t/ 2m (t), it is easy to see that 

(2.16) U 2m _2j(tfc) = U 2m _2j(l?2fc,2m) = —U2j-l{jl2k,2m) = — U 2 j_i(tfc). 

Indeed, since sin(2m + l)p 2 k, 2 m = 0 and cos(2 m + l)i) 2 k, 2 m = 1, it follows from 
the addition formula that sin(2m — 2 j + 1 )i] 2 k, 2 m = — sm 2 jr) 2 k, 2 m, from which the 
equation follows. 

Since U 2m (fi) = ... = Z7 2m (t m ) = 0, we evidently have 
detEj(t) = U 2 m(to) detEj(t), 

where E,- (t) is a submatrix of E j (t) with the column Xj (t 0 ) and the row containing 
U : 2 m{tj) of Ej (t) removed. Using the equation 112.161) . we can replace all rows of 
Ej (t) that contain even Chebyshev polynomials by rows that contain odd Cheby¬ 
shev polynomials. More precisely, we replace the rows (t/ 2 ,;(U),..., f7 2 i(f m )) for 
i = j,j + 1, • • • ,m — 1 by rows (f7 2m _ 2 j_i(ti),..., U 2m _ 2i _i(t m )), respectively. The 
new matrix has all elements given in terms of the Chebyshev polynomials of the 
odd degree; hence, 

detSj(t) = edet(f7 2i _i(4))™’ 1 ™ =1 

where e = ±1. Assuming now that this determinant is zero. Then there exists a 
non-zero polynomial, Q, of the form 

m 

Q(x) = 'Y' J b i U 2 i-i(x), 

i= 1 

which vanishes at t\, ..., t m . Since Q is odd, it vanishes also at — 1 \,..., —t m . Ob¬ 
serve that the sets (— 1 \,..., —t m ) and (ti,..., t m ) do not overlap. We conclude that 
the polynomial Q of degree 2m— 1 vanishes at 2m distinct points and, consequently, 
it vanishes identically, a contradiction. □ 

A similar theorem holds for n = 2m — 1, which we state below. 

Theorem2.12. Lettj = r/ 2 j, 2 m for j = 1,2,... ,m. Then the matrices Ei(t), E 2 (t), 
. ..,S TO (t) in Theorem 2.1 are all nonsingular. Consequently, the polynomial P £ 
H 2m _i is uniquely determined by the set (im of its Radon projection. 
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The proof is similar to that of Theorem 12.1 H and it uses the equation 

U2m-2j{V2k-l,2m) = Lgj-1 i?l2k- 1,2m), 

which can be verified using elementary trigonometric identities. 

We have conducted some numerical tests for identifying other sets of t k for 
which the matrices 3i(t),..., 3 m _i(t) are nonsingular, so that the reconstruction 
of a polynomial from the set of its Radon projections is unique. For n < 20, 

it turns out that both the equidistant points in (0,1) and the Chebyshev points in 
(0,1) work out. See, however, the discussion at the end of the next section. 


3. Reconstruction of polynomials 


In this section we set up the algorithm that can be used to compute P n from the 
Radon projections. We only consider the case n = 2to. Let ~{j. k be the data 

(3.1) 7 j,k = (/; tk)/\J 1 - if, 0 < j < 2m, 0 <k<m, 

where 4>j, m are given as in 12.1311 and tfc are distinct numbers in [0,1), chosen in 
advance, such that the linear systems of equations (IQ) and iTHl) below have unique 
solutions for all j. 

Recall that the Lagrange interpolation by trigonometric polynomials is used in 
the proof of Theorem l2.(il The Lagrange interpolation based on the points in (ETT31) 
is given explicitly by (see, for example, mi) 


2m 

L n f{4>) = 


3=0 




sin(m + \){<t> - </>y m ) 
(2m + 1) sin \{(j) - ' 


where we assume that / is the function being interpolated; that is, L n f(4>j trn ) = 
for 0 < j < 2m. Using the well-known formula 

sin(m + i )6 

1 + 2 cos (/) + ... + 2 cos m</> =-r-, 

sinf 

we can write L n f in the standard form of a trigonometric polynomial, 

L nf{4>) = m 0 (f) + 2^2 (mf (f) cos j(/) + mj(f) sin j(f>) 

3 =i 


mf(f) = fi&l,™) COS 

m + 1 1=0 
^ 2m 

m S j (/) = 2-n+i fi^rn) sm 

For each k, 0 < k < m, we will use the above formulas with ) = "fj,k and we 

shall write rn^ k = mf(f) and k = m^{f) for this particular /. This will allow 
us to determine the values Ao(t k ), Aj(t k ) + A 2 m -j+i(tk) and Bj(t k ) — B 2 m -j+i{tk ) 
for 1 < j < m. 

The next step is to fix j and use the values of Ao(ffc), or Aj(tk) + A 2 m -j+i(tk), 
or Bj(tk) — B 2 m-j+i(tk) for k = 0,1 ,..., to to determine the coefficients of Aj and 
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Bj. For this purpose we consider the following linear systems of equations: For 
1 < j < rn/2 , 

m m 

(3-2) E dj,2lU2l{tk) + E dm—j-\-l,2l U 2 i-\{tk) = m 2 j,k, 0<k<m, 

1=3 l=m—j+1 

and for 1 < j < (to + l)/ 2 , 

m m 

(3-3) E dj,2l-lU2l-l(tk) + E dm—U 2 i(tk) = m 2 j-i,k, 0 < fc < TO. 

1=3 l=m—j +1 

The quantities TOj.fc will be specified later. It is easy to see that the coefficient 
matrix of these systems of equations are S j in Theorem 12 .(il There are total to . 
system of linear equations, each of size (to. + 1) x (to + 1). Solving these equations 
with proper nij : k will determine the coefficients of Aj and Bj. 

The last step is to use (E3 and (l 2 ~ 8 l) to get Cj t k, which are the coefficients of P n 
in ( 12 . 211 . 

Algorithm: 


Step 1. For 1 < l < to and 0 < k < to, compute 


(3.4) 


2 m 

V 


ms l 


and 


™?,k 


i 

2 TO + 1 


2 m 

j,k sin Wj. 

3=0 


Step 2 . Solve the systems of equations (14.211 and (14.411 to get the coefficients (ij.k 
and bkj in (ESI) : 

Case 1. Solve the to systems of ESJ and (TH 1 ) for nij } k = rn^ k to get 

ayfc := dyfc, 1 < j < to 0 < k < 2 to. 

Case 2. Solve the to systems of 14.211 and 14.411 for nrij^ = rrij k to get 

. . to + 2 

bj, 2 i ■= djji, 1 < 2 < m/2, and b ji2 i ~-d jt 2 i, —^— < 3 <m, 

TO + 1 TO + 1 

0j,21-1 ■= dj ,21-1, 1<J<—2—’ and b J,21-1'■=-dj,21-1, —2— 

where 1 < l < to. 

Step 3. Substitute the output aj,k and bj, k in Step 2 into E3 and ESI to get 
Cj./j. Then the polynomial P is given by ESI) . 


The output of this algorithm is the polynomial P that satisfies 03J. The main 
computation appears to be the Step 2. Initial numerical experiments indicate that 
the linear systems of equations ESI and ES1> are ill-conditioned. However, among 
the set of points that we tested, the largest condition number of the matrices 
corresponding to the equidistant points t k = (k + l)/(m + 2 ), 0 < k < to, of ( 0 , 1 ) is 
smaller than that corresponding to the Chebyshev points tk = cos(fc + 1 ) 7 t/( 2 to+ 4 ), 
0 < k < to or the points in Theorem mm 

As pointed out by one referee, it is perhaps not surprising that the matrices of the 
linear systems ESI and ESI) are ill-conditioned. Solving these linear systems means 
inverting the Radon transform, while it is known that the Radon reconstruction 
problem is ill posed. Thus, for to large, the algorithm may not be useful for practical 
computation. On the other hand, our result shows that it is possible to determine 
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a polynomial of degree n uniquely from a set of (n + l)(n + 2)/2 Radon data that 
concurs with the parallel geometry. This appears to be of independent interest. 
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